Do correlations create an energy gap in electronic bilayers? 
Critical analysis of different approaches 
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This paper investigates the effect of correlations in electronic bilayers on the longitudinal collec- 
tive mode structure. We employ the dielectric permeability constructed by means of the classical 
theory of moments. It is shown that the neglection of damping processes overestimates the role of 
correlations. We conclude that the correct account of damping processes leads to an absence of an 
energy gap. 
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Layered electronic systems have been of intense interest recently. Such systems can now be routinely fabricated 
in semiconductors with well-controlled system parameters. We focus in this paper on the analysis of the collective 
i— - n l modes of electronic bilayers. These systems consist of two quasi-two-dimensional (2D) layers of electron liquids with 
electron density n s , separated by a distance d comparable to the interparticle distance a — {ris/ir) 1 / 2 within the 
layers. Quasi-two-dimensional here means that the electrons have quantized energy levels along one dimension, but 
are free to move in two dimensions. A well known example of 2D electron system are electrons confined in the 
vicinity of a junction between a semiconductor and isolators [in a MOSFET structure] or between layers of different 
semiconductors [in heterojunctions]. In these systems electrons are confined near the surface by an electrostatic field; 
in electronic bilayer systems two electron layers are separated in a double quantum well. Another class of layered 
systems, electronic superlattices, consisting of a large number of identical electronic layers, show a similar behavior 
\ and are beyond the scope of this paper. High r s (r s — a/as] clb — s s h/e 2 m* being the effective Bohr radius) values 
£^ ■ are now available in 2D layers jlj and the technique should be available to fabricate relativley high r s bilayers also. 
One expects that at high enough r s values (r s > 20 0) a bilayer crystallizes into a Wigner lattice. Here, we restrict 
our considerations to the strongly coupled liquid phase bilayer systems. There are some theoretical investigations in 
this strongly- coupled liquid phase regime, both on the static and the dynamic level ]3|-^10||. 
t-H , One of the issues under consideration is the analysis of collective excitations in bilayers. First we recall the Random 
phase approximation (RPA) result applicable in the weak coupling regime r, C 1. In the RPA there are two 
longitudinal modes: an in-phase and an out-of-phase mode. In the first mode the two layers oscillates in phase and the 
resulting dispersion relation is similar to that of an isolated 2D layer; for k — > the eigenmode behaves asw~ Vk. In 
1^ \ the out-of-phase mode the oscillation phase of the two layers differs by tt and we deal with an acoustic mode u> ~ k as 
k — > 0. However, in the strong coupling regime r s 3> 1 the collective modes structure is strongly affected by particle 
Q\ ' correlations. The effects of correlation beyond RPA can be recasted into the local field correction. Early approaches 
focused on intralayer correlations but ignored or suppressed interlayer correlations 0. There have been two basic 
lines which also take into account the interlayer correlations. One line uses a low frequency (static) approximation for 
the local field correction. This can be done by applying the Singwi-Tosi-Land-Sjolander (STLS) approximation |l^] 
to the bilayer problem The other line uses a high-frequency local field correction. The analytical most simple 

expressions can be obtained using the quasilocalized charge (QLC) approximation Related studies to the bilayer 
system have been done in Refs. ||[|. Examiming the effect of the particle correlations on the collective excitations 
structure the two methods have arrived at different results. In particular, the QLC method predicts the occurence of 
. a finite energy gap (u> > for k = 0); however, no such energy gap appears in the calculations based on the STLS 
approximation. In Ref. [Q the formal reasons leading to the different results are clarified. It has been shown, that in 
. i-H ! the STLS formalism the local field correction vanishes as k — > and as a result we arrive (qualitatively) at the RPA 
dispersion expression without an energy gap. On the other hand, the high-frequency QLC local field correction docs 
not vanish at k = and this leads to the occurence of the energy gap. This paper is adressed to the question whether 
. - - i interlayer correlations create an energy gap in bilayer systems. 

The bilayer, consisting of two 2D electron layers of area A embedded in a neutralizing background, and separated 
by distance d can be mapped onto a two-component two-dimensional system |6|. In the following we neglect the 
layer thickness, the influence of metallic electrodes and the difference between the dielectric constants of different 
semiconductors. Then the corresponding interaction potentials are 

<pu(k) = <f22{k) 



k ' 

<p 12 (k)= 2 -^e- kd . (1) 
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The two component system may be described by a matrix formalism in species space p^qjlq ] . However, in the present 
case of two identical electron layers the corresponding matrices diagonalizes and it is much easier to investigate the 
scalar dielectric functions for the in-phase and out-of-phase motions, £,„ and e out respectively, 

e^ 1 (k,w) = 1 + <PaXa(k,a;) , a = m,out. (2) 

where ¥>i n ,out = 27re 2 (l±e~ M )/fc are the interaction potentials for the in- and out-of-phase motions, respectively, and 

X a (k,u/) = (HA)- 1 f dv [ dt < [n Q (r,i),n Q (0,0)] > ^ T+iut , (3) 



are the density response functions. Here ni n ,out(r> t) = n>i(r, t) ± n 2 (r, t), n-i(r, t) (i = 1, 2) being the electron number 
density operator of the «th layer, [A, B] is the commutator of operators A and B, and < A > is the equilibrium 
average of operator A. 

The calculations based on the QLC formalism jl3|,^| lead to the following expression for the bilayer dielectric 
permeabilities, 

x w\{k){l + e- kd ) 

£ in ( k > w ) = 1 + 73 



2 - wf(k)(l + e-W) - (D n + D 12 ) ' 
-m s 1 w 2 (k)(l - er kd ) 

cj^ — wf — e Ka ) — (Dn — D12) 

where Wi(k) = {2ire 2 n s k / 'm) 1 / 2 is the 2D plasma frequency of an isolated electron layer with surface number density 
n s . The positive magnitudes Dn and _D 12 take into account the intra- and interlayer Coulomb correlations between 
the electrons, and are expressable via the Fourier transform of the pair correlation functions hn(k) and hi2(k), so 
that HQ 

D 11 (k) = -w 2 1 (k)J2 { ^- {M|k-q|)-Mff)-e- ,d M<*)} - 

q 

D 12 (k) = w 2 (k) ]T h 12 (\k - q|) . (5) 

q 

The poles of the inverse dielectric permeabilities determine the eigenfrequencies in the bilayer system. We find in 
accordance with |]|| 

^ 2 n (fc) = wl(k){l + e- fed ) + (Dn + Dn) , 
ul ut (k) = w 2 (fc)(l - e- fed ) + (D n - D 12 ) . (6) 

These expressions were analyzed in Refs. It was found that the in-phase modes are not qualitatively different 

from the corresponding modes in the isolated 2D layer (with double density). In particular, for k — > the typical 
2D soft plasmon mode behavior uj ~ \fk can be observed. On the contrary, from Eq.(^) one concludes that the 
out-of-phase modes develop an energy gap at k = |6|,[|[l4| : 

w 2 (0) = _^5/ dqq 2 e -«<l hl2 (q) . (7) 



II 



However, this result looks quite strange. One expects that in the longwavelength limiting case k — > the wave 
does not feel the separation between the two layers and we should deal with a double density single 2D layer. For 
the single layer, however, an out-of-phase mode is unknown. Moreover, as can be seen from Eq.(Q) the energy gap for 
the out-of-phase mode remains even for the case of a vanishing layer separation d — 0, i.e. for the case of a single 2D 
layer. Due to the translational invariance of the double density single 2D layer only the plasmon mode oj 2 (k) — 2wf(k) 
should be observed as k — > 0. One concludes therefore that there should not be any energy gap in the bilayer system. 
Summarizing the above discussion, we have found that Eq.(^) predicts an energy gap which should not exist. 

To clarify this discrepancy we go back to Eqs.(H). These equations describe the dielectric permeabilities of systems 
of oscillators with eigenfrequency tof a (k) (or uj 2 ut (k)) and with the oscillator strength f in = l + e~ kd (or / out = 1 — e~ kd , 
respectively). We focus now on the out-of-phase mode and consider the long- wavelength limiting case k — > 0. The 
eigenfrequency of the out-of-phase mode reduces to the energy gap value (Eq.(^|)) but the oscillator strength tends to 
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zero as / ou t = kd. This is just the solution of our problem. Within the QLC approach there is indeed a finite energy 
gap for the out-of-phase mode, but this mode does not develop as k — > since the mode oscillator strength vanishes. 
In a quantummechanical language this means that the transition between two energy levels separated by an energy 
gap is forbidden; in a classical language one would say that the number of oscillators having the eigenfrequency w 2 (0) 
(Eq.(g)) tends to zero as k — ► 0. 

Our discussion results in the following picture within the QLC approach: as k — ► the eigenfrequency of the 
out-of-phase mode tends to the energy gap value but the mode signal becomes weaker and weaker. Now a deficancy 
of the QLC algorithm comes into play. The QLC method is not able to describe damping processes and therefore 
the calculations based on QLC approximation fail to provide information on the damping of the collective modes. 
However, with decreasing mode signal the noise to signal ratio increases and it is necessary to take into account 
damping processes. 

In what follows we employ the classical method of moments to determine the dispersion law of the bilayer collective 
modes. To do so we use an interpolation formula, satisfying all known exact relations and sum rules. This allows 
us to include damping processes into our calculation, at least formally. In Ref. |l7| an expression for the dielectric 
permeability of a single 2D electron layer satisfying all known sum rules is constructed. Generalizing the results 
of p7| to the present case one expresses the dielectric permeability for the out-of-phase motions via the frequency 
moments Mi = — — u* Im£ ( 7 u 1 t (k, oj) du> (i = —1,1,3) and the function g out = g ut(k, z) being analytic in the 
upper half-plane Im z > , z = lu + i rj and having a positive imaginary part there, and such that [q on t{k, z)/z] — > 
for z —* oo, 

-x f , s wl(k)(l - e- kd ) (o; + gou t(k,a;)) 

e outl k > W J - 1 H 7 \ 7 — Y ■ (°) 

loLo* - U7?(fc)(l - e- kd ) - (D u - D 12 )j + q out (to 2 - w\{k){l - e - kd ) [l - e^ t (k)] 1 J 

Eq.(||) interpolates between the exact high-frequency behavior (given by the first and third frequency moment sum 
rules and the low- frequency behavior (determined by e ou t(k) - the static dielectric permeability for the out-of-phase 
motion). For simplicity, in Eq.(Q) we have omitted the kinetic energy term into the third frequency moment since 
we are interested in the strong coupling limit only. An expression similar to Eq.(||) can be obtained for the in-phase 
motion. Notice that the QLC expression Eq.(|^) can be obtained by putting q out = 0. The QLC expression also satisfies 
the first and third frequency moment sum rules but does not reproduce the low-frequency behavior of the dielectric 
permeability. Via the inclusion of a finite (complex) function q out one is able to describe the low-frequency behavior 
in the correct manner. In addition we have now the formal possibility to include damping into our considerations. 
In the liquid phase of a strongly coupled bilayer the damping is mainly given by the diffusion of the quasilocalized 
particles site positions. A nonzero diffusion constant for the site positions migration results in a nonzero conductivity 
for the out- and in-phase motions. We mention here that the 2D static "conductivity" cr ou t is not a real conductivity 
since it describes the out-of-phase flow. Notice that the existence of a nonzero static out-of-phase "conductivity " is 
guaranteed only if q ou t(k — > 0,0) = ih, with h = 4ira ou t(Dn — D12) / diw\(k)(l — e~ kd ) ~ k~ 2 (di being an effective 
thickness of one layer). We have no phenomenological way to choose the exact q(k, uS). However, the most natural 
way to generalize the QLC result is to put ? ut(k,w) equal to its static value in the long- wavelength limit q out = ih. 
Within this approximation the high frequency pole of the out-of-phase permeability is shifted into the lower complex 
half plane. Instead of Eq.(@) we have now: 



h / h 2 

u J 2 out (k^0) = - t -± ] l~- + (D 11 -D 12 ) , (9) 

and the mode corresponding to the energy gap in the QLC approximation becomes overdamped. 

In the previous discussion we have shown that the interlayer correlations lead to a divergence of the third to the 
squared first frequency moment ratio M^/Mf in the long-wavelength limit. Using the fact that q ou t ~ M3/M 2 the 
expansion of Eq.(||) in terms of Mf /M3 provides as k — > the following expression: 

e- u \(k - 0, w) = 1 + w\{k){l-e- kd ) + a _ 

lu 2 + Luq'(k, u) - wl(k)(l - e-W) [1 _ £ "i (fc)] 

Here q'(k,Lu) = —(Du — Di2)/Q'out(k, u>) is a function which has the same properties as <7 ou t(k, uj) has. Notice that 
Eq.(|l0|) contains the information on the static dielectric permeabilty e~ u 1 t (fc), which is given by a compressibility sum 
rule as k — > 0. The high- frequency part is given only by the first frequency moment, whereas the "diverging" third 
frequency moment is absent. Eq.(|lO|) corresponds therefore to a STLS-like approximation for the out-of-phase motion 
dielectric permeability. We have again no phenomenological way to choose the exact q'(k.,uj). However, from Eq.(|i"o|) 
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it can be seen that v = — iq'(k, ui) can be regarded as an effective collision frequency in a Drude-Lorentz like theory. 
Neglecting the frequency dependence of this collision frequency for the out-of-phase motion one defines the collision 
frequency by the static "conductivity" of the out-of-phase motion |l7j er ut, v(k) = w\{k)(l — e~ kd )di/4ira out . We 
will not perform a detailed analysis of this quantity. We mention here only that in a case of a finite out-of-phase 
"conductivity" we obtain that v{k) behaves as k 2 for k — ► 0. Then from Eq. (|To|) we obtain the following dispersion 
for the out-of-phase motion in the long-wavelength limiting case, 

^ {k) = ^Ae +0 (k 2 ). (11) 

TO 

This is just the result of the STLS calculations |l4|]. Notice that Eq.(|TT|) corresponds to the low- frequency pole of 
Eq.(||). Thus we have found that the expression Eq.(||) calculated from the theory of moments contains the QLC 
expression and a STLS like expression as limiting cases. Further we have shown, that in the long wavelength limiting 
case Eq.(^) converts into a STLS like expression and does not predict the existence of an energy gap in the bilayer. 
Due to the absence of an energy gap in the present approach we observe a soft transition from the bilayer to the 
double density single layer with layer separation d — > 0. However, within the qualitative analysis of this paper we are 
not able to establish the value of wavevector k at which the transition to the STLS like permeability occurs. However, 
some qualitative statements can be made. The "transition point" to the STLS regime is mainly given by the function 
Qout ~ v~ 1 (k). With increasing coupling strength the diffusion constant decreases and so docs v~ 1 (k). We also expect 
that the STLS "transition point" is shifted to lower k values and it might be possible to observe an "energy gap" at 
not to low k values. 

In this Brief Report the dispersion law for the out-of-phase mode in an electronic bilayer is analyzed. The analysis 
is based on the expression for the dielectric permeability of the out-of-phase motion obtained from the classical 
theory of moments without using perturbation theory. We have shown that this expression contains the permeability 
expressions obtained by other approaches as the QLC or the STLS approximation as limiting cases. The analysis of 
our expression has confirmed the STLS predictions of an absence of an energy gap in the bilayer system. Finally, 
we invite the experimentalists to verify the different predictions of the various theoretical approaches concerning the 
existence of an energy gap in a bilayer. This would exhibit not only a check for this specific system but also a general 
verification of the various theoretical approaches applicable to other plasma systems. Another tool to investigate the 
collective mode structure in bilayers might be numerical simulations. However, due to the finite number of particles in 
the simulations the available k values have a lower bound. Therefore one remains with the difficult task to extrapolate 
the behavior at k — from the behavior at relative high k values. 
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